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We introduce a procedure to measure the density matrix of a material system. The density matrix 
is addressed locaUy in this scheme by applying a sequence of delayed light pulses. The procedure is 
based on the stimulated Raman adiabatic passage (STIRAP) technique. It is shown that a series of 
population measurements on the target state of the population transfer process yields unambiguous 
information about the populations and coherences of the addressed states, which therefore can be 
determined. 
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Active manipulation of the quantum state of different 
microscopic systems increases the need for measurement 
schemes which verify the reliability and efficiency of the 
engineering procedure. Quantum state reconstruction 
methods have been proposed in several systems, includ- 
ing light field and matter systems such as the vibrational 
state of molecules, trapped atom motion, Bose-Einstein 
condensates, atomic matter waves, electron motion etc., 
for recent reviews see Refs. 

In this work we introduce a procedure to measure the 
density matrix of a material system. The density matrix 
is addressed locally in our scheme. More precisely, the 
measurement process yields a part of the total density 
matrix, namely the elements 



(1) 



where the indices m and n refer to some steady states of 
the material system. The procedure is based on popu- 
lation measurements, so one needs an ensemble of iden- 
tically prepared systems to obtain the required density 
matrix elements. In practice this means that the mea- 
surement is performed for example on an atomic or a 
molecular beam. A population measurement inevitably 
entails a reduction of the quantum state, therefore the 
state of the system is destroyed in the process. The 
complete density matrix can be determined by repeated 
application of the measurement procedure for different 
pairs of states {?7i',n'}. Recently, some new proposals 
have been published to achieve a similar goal [j. Hi • We 
will compare our method with those works at the end of 
this paper. 

Our measurement procedure is based on the stimulated 
Raman adiabatic passage (STIRAP) process Its im- 
plementation imposes the following requirements on the 
system: the linkage displayed in Fig.^must be realizable. 
We assume that the states labeled by |2) . . . , |A^)} 
are populated, the others are empty. There are three 
light pulses fri{i), fa (^)} which couple the states 

{|m), |a)} via an excited state |e). The other popu- 



lated states must remain unaffected by these light fields. 
This could be achieved by exploiting selection rules or se- 
lection based on resonance frequencies. The frequencies 
of the light fields may be detuned by A from the transi- 
tion frequencies Lo^i = (i?e — Ei)/h, where Ei is the energy 
of the state |i) G {|a), . . . , |A^}}, but the three-photon 
resonance condition must be fulfilled. 

Initially, the states \a) and je) are empty. The Master 
equation which describes the measurement process reads 



dg 
dt' 



[H{t),Q]~^{\e){e\g + g\e){e\) 



r, 



^{\a){a\g + g\a}{a\), 



(2) 



where g denotes the density operator of the system. The 
constants and Fa stand for the decay rates from the 
states |e) and |a), respectively. The Hamiltonian H{t) 
is given in the interaction picture and the rotating-wave 
approximation (RWA) as 

H{t) = hA\e){e\ + ^ J2 mm){e\ + h.c), (3) 



where the pulsed Rabi frequency ^i{t) derives from the 
field £i{t). The Rabi frequencies il,n{t) and iln{t) are 
taken to have the same envelopes: 



^m{t) = ^p{t) COS a, ^n{t) = ^p{t) sina e 



t(3 



(4) 



where a and /3 are fixed angles which define the relative 
amplitudes and phase of the two pulses, respectively. The 
pulses TO, n and a are delayed with respect to each other; 
however, for an efficient STIRAP process they must have 
a significant overlap. In the following Qp(t) and Qa{t) will 
be taken real. 

Let us neglect the dissipative terms in the Master equa- 
tion Then we have a purely unitary evolution. Now, 
our aim is to determine the time evolution operator U{t) 
which governs the time development of the density oper- 
ator in the absence of dissipation. It will help us under- 
stand the essence of the measurement procedure, and we 
later return to the discussion of dissipations. 

The measurement procedure consists of performing a 
STIRAP process in which a part of the population from 
the states |m) and \n) is transferred to the state \a). The 
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pulses m,n, which act on the states {|to), |n)}, play the 
role of the pump pulses, whereas the pulse a corresponds 
to the Stokes pulse. It is assumed that they arrive in 
the counterintuitive time order. The pump pulses m,n 
define a coupled state (or bright state) 

|C) =cos a |m) + sina e'^'ln), (5) 

and an orthogonal decoupled state (or dark state) 

|Z)) = — sin a |m) + cos a e*'^|ri), (6) 

see Ref . . The Hamiltonian in Eq. (PJ can be expressed 
in terms of the states |C) and \D) to obtain 

H{t) = hA\e){e\ + ^{np{t)\C){e\ + na{t)\a){e\+h.c.). (7) 

Obviously, the Hamiltonian docs not act on the decou- 
pled state \D). From this formulation it can clearly be 
seen that we have in fact an ordinary STIRAP process 
defined on a three-state system. The dark state \4>o{t)) 
of the Hamiltonian in Eq. |2J is given by 

l^Polt))^ cos e\C) -sin 9\a), (8) 

and the two bright states are 

|'0+(t)) = sin(/?sin6'|C) -|- sin(pcos0|a) + cos(y5|e) 

'(/;_(t)) = cos(^sin0|C) -I- cos(^cos0|a) — sin(^|e), 
where the angles 9 and ip are defined as 



(9) 



tan( 



tan 2(p- 



n{t) 



(10) 



with n{t) = ^^l{t) + i^l{t), see Ref. 5]. In the adiabatic 

limit, a simple calculation shows that the unitary time 
evolution operator U (t) reads 

U{t)= \D){D\ + \Mt)){C\ 



+ cxp ( / e+{t')dt' ] |^+(t))(^+(-oo)| 



+ exp ^~-J s^it')dt'j |^_(i))(V'-("0<)|P,) 

where e+{t) — ^^{t) cot tp and e-(t) — —^fl{t) tan tp are 
the two nonzero eigenenergies of the Hamiltonian Eq. ||7J). 
The physical interpretation of this result is quite simple: 
The decoupled state \D) remains untouched throughout 
the transfer process. The coupled state \C) follows adia- 
batically the dark state \tpo{t)) of the Hamiltonian H{t) 
[Eq. O]. The bright states \i^±{t)) also evolve adiabat- 
ically, however, they acquire a phase shift because they 
belong to nonzero eigenenergies e± (t) . 

In the adiabatic limit, the bright states \ip±{t)) are 
not populated throughout the whole time, provided that 
initially they were not populated. Thus, we have found 



that the final state of the system, after the pulses have 
passed becomes 

lim U{t)g,U^it) = \D){D\g,\D){D\ + \a){C\g,\C){a\ 

- \D){D\g,\C){a\~\a){C\g,\D)m) 

Let us suppose that we are able to measure the pop- 
ulation Pa on the state \a). Then, the result can be ex- 
pressed in terms of the density matrix elements {gi)mrm 

{gi^nn^ {gi^mn aS 



Pa = {C\gi\C) = cos^a(ft)„„i H-sin^Q;(ft)r 
-I- sin2Q;Re{(ft)„„e*^}, 



(13) 



where we have made use of the definition Eq. ||SJ). The 
outcome of the measurement is phase sensitive: It de- 
pends on the relative phase (3 between the fields m and 
n (cf. Eq. I0J). Observe, that the coherence {gi)mn also 
appears, so one expects that by accomplishing an appro- 
priate series of measurements one can obtain not only the 
diagonal elements of the density matrix, but the coher- 
ences as well. Indeed, the three required density matrix 
elements can be determined in four steps: In the first 
and second measurements, one of the pump field n or 
TO is switched off which corresponds to taking a = or 
a = 7r/2. In this way, we get the diagonal density ma- 
trix elements {gi)mm and {gi)nn- In the third and fourth 
steps, both of the pump fields to and n are present so that 
the value of a is chosen as 0<a<7r/2. One measure- 
ment is performed with the choice f3 — and an other 
with (3 = —Tr/2. It is easy to show, that by combining 
the results of the last two measurements with that of the 
first and second, one can unambiguously determine the 
needed (ft)mn- 

There is one question left, which we have to answer: 
How can one measure the population Pq? Here we return 
back to the analysis of the dissipation mechanisms which 
we have postponed so far. Let us consider again the 
Master equation Q . If the usual STIRAP conditions are 
met during the measurement procedure, then the decay 
from the excited state |e) has got negligible effect, since 
this state is minimally populated 7, 8]. The situation is 
quite different in the case of the decay from the target 
state I a): This state is deeply involved in the transfer 
process. The dissipation may influence signiflcantly the 
population transfer to the target state and as a result, 
it deteriorates the fidelity of the measurement. However, 
it can also be utilized to monitor the population on the 
target state \a). 

It is plausible to assume, that the decay rate from the 
state I a) to the states |to) and |n) is negligibly small, since 
they have the same parity, because both of the transitions 
|m) <-> |e) and |e) <-*■ \a) are dipole allowed; the argument 
is the same in the case of the state \n). The decay from 
the state |a) to the state |e) is undesirable. It can be 
avoided if the energy of the state \a) is lower than that of 
the state |e) or if there are faster decay channels. In the 
following we assume that the decay between these states 
is negligible. 
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Now we present a simplified treatment of the dissipa- 
tion from the target state \a) by assuming that the ini- 
tial state of the system can be given by a state vector 
B = [Be, Be.Bay in the basis {|C), |e), \a)). Note that 
in this definition we have included only those states of 
the system which participate in the population transfer 
process. This means that the norm of B may be smaller 
than unity. The Master equation Q is replaced by a 
Schrodinger equation with a non-Hermitian Hamiltonian. 



J 



In the adiabatic basis Eqs. 
equation reads 



and lO, the Schrodinger 



■H'{t)A{t), 



(14) 



where the effective non-Hermitian Hamiltonian H'(t) is 
given by 



H'it) = 



cot (f — iFa cos^ 9 sin^ ip i (^^Ta sin(26') + tij sni ip 
i (^Ba sm{2e) - 0) sin (p -iPa sin^ 9 



'ijBa cos^ 9 sin(2(/5) + i(p 
(ir,sin(20)- 



cos ip 



i^Fa cos^ 9 sv[i{2pi) ~ ip i ( ir'a sin(26') + ^) cosi/j ^VHanip — iFa cos^ 9 cos^ p 



(15) 



r 



Decay from the excited state |e) is neglected in order to 
obtain more simple expressions. This approximation is 
justified by the small involvement of the excited state and 
by assuming small decay rate Ff.. The state vector B{t) 
is transformed by the orthogonal rotation 



0{t) = 



according to 



sin Lp sin 9 cos 9 cos p sin 9 

cos p Q — sin p 
sin ip cos 9 — sin 9 cos p> cos 9 



A{t)=0-^{t)B{t). 



(16) 



(17) 



In Eq. (|15|l one can see that the decay from the target 
state I a) affects the time evolution both of the dark state 
and the bright states. More precisely, terms proportional 
to Fa appear not only on the diagonal but also on the ofF- 
diagonal part of the Hamiltonian. On one hand, these 
terms lead to a decrease of the probability amplitudes 
for all of the three adiabatic states; on the other hand, 
they contribute to nonadiabatic couplings which mix the 
adiabatic states among themselves. If r^max ^ A and 
the process is adiabatically slow, then the diagonal ele- 
ments dominate over the nondiagonal ones in Eq. 
Then the states \il)+{t)) and |V-'-(^)} can be adiabatically 
eliminated by setting A+(i) = A_ (i) = in Eq. and 
solving the resulting algebraic equations for A+{t) and 
A^{t). Thus, by inserting these solutions in the differen- 
tial equation for Ao(t) we find the final population Pa{Fa) 
on the state \a) to be 



Pa{ra) 



Poexp 



-2Fa 



A2g2 cos^g + (»V4 + y)^)^sin^ 

A2r2cOs4 6'-K (f^V4-Fv?2)2 



(18) 
-dt 



since Pa{Fa) = |Ao(oo)|2. The symbol Pa denotes the 
initial population on the state |C) which is given by 
Eq. (|13|l . In the absence of decay Fa = we recover 
the result derived earlier. For a significantly large decay 
rate Fa, the target state \a) is emptied by the end of the 
population transfer process. If the decay is accompanied 
by emission of photons, then the signal will be propor- 
tional to the population Pa- 
in general, the initial state of the system cannot be 
described by a state vector. In this case, our consid- 
erations in the preceding paragraphs cannot be applied 
straightforwardly. However, some results are still valid: 
In the Master equation ^ the projection of the den- 
sity operator g on the decoupled state \D) is a conserved 
quantity during the transfer process. This implies that 
the maximum amount of population which can be trans- 
ferred to the target state \a) is given by Pa in Eq. H13|l . 
Therefore, we have performed a numerical simulation, 
whose main purpose is to verify whether this portion 
is indeed transferred. The simulation may also allow 
us to study the adiabaticity of the process: the excited 
state must be involved only minimally in the population 
transfer. As a result we have verified numerically, for 
a wide range of Fa and for several initial density oper- 
ators gi, that the population (C|p/|C) is zero and the 
population on the excited state |e) is very small through- 
out the whole time. The quantity {C\Qf\C) is obtained 
from the simulation in the presence of decay from the 
state I a). The parameters of the simulation have been 
chosen as follows: The atomic system is described by 
{A,Fe,Fa) = (0.3,0.1,0.0 - 3.0); the pulses have been 
Gaussian with parameters (fimax, ''') = (6.0,2.0,3.2), 
where T denotes the half-width of the pulses, and t is 
the pulse delay. Time and frequency are measured in 
arbitrary units. 

In this way we have shown that combining a STIRAP 
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process with a population measurement on the target 
state of the population transfer, yields the expectation 
value of the operator 

V = cos^ a I to) (m I + sin'^ a I n)(n I 

+ sina cosa(e*'^|TO)(n| + e"'''^|n)(TO|). (19) 

Performing a series of measurements of this kind by vary- 
ing the angles a and /?, the density matrix elements 
(ft)mm, and (ft)mn Can be unambiguously ob- 

tained. In Refs. a measurement procedure has been 
proposed where the population from the ground state en- 
semble is transferred to an excited state. This scheme 
seems more sensitive to the decay from the excited state 
because it gets significantly populated during the pro- 
cess. Therefore, the incoherent decay back to the ground 
states may interrupt the coherent evolution which is re- 
quired in this measurement procedure as well. 

In conclusion we have worked out a measurement pro- 
cedure based on stimulated Raman adiabatic passage to 
obtain the the density matrix elements of a material sys- 



tem. The scheme is highly reliable due to the robustness 
of the STIRAP process. In the beginning of the measure- 
ment, a part of the density operator is expressed using 
two orthogonal states: the coupled and the decoupled 
states, which are defined by the two pump pulses. In this 
way a 2 X 2 block of the density matrix is addressed. The 
coupled part of the density operator is then transferred 
to an auxiliary state. We have shown that by measuring 
the population on the auxiliary state for some configura- 
tions of the pump pulses, the 2x2 block of the density 
matrix can be determined. 

The measurement procedure presented here can be 
used efficiently when it is not necessary to obtain the 
complete density operator, but only few elements are 
needed. We propose the implementation of this scheme 
in several microscopic quantum systems, e.g. electronic 
states of atoms, vibronic states of molecules etc.. 
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FIG. 1: Interaction scheme for the measurement of the den- 
sity matrix. The states labeled by arabic numbers are occu- 
pied, the others are empty initially. The three laser pulses 
£i{t) couple two of the populated states |m) and |n) with the 
auxiliary state \a) through the excited state |e). A common 
detuning A is allowed but the three photon resonance is re- 
quired. The wavy arrows indicate decay from the excited and 
auxiliary states. 



